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Abstract 

We perform a chiral one-loop calculation of the unitarity corrections to the processes K + — ► 
7r + 7/ + / _ up to C(E 6 ), taking into account 7r + 7r~ intermediate states. Branching ratios and differ- 
ential branching ratios are computed and presented to demonstrate the importance of the above 
corrections. 

I. Introduction 

The investigation of radiative rare kaon decays has taught us that they form a complex of interrelated 
processes which share some common features. Experimental and theoretical results on any of these 
reactions are useful in the analysis of all of them. Since they can be analyzed using chiral perturbation 
theory (ChPT) HI], the experimental exploration of the entire complex provides stringent checks on 
this theoretical method. Recently, radiative kaon decays have attracted considerable attention from the 
theoretical |2], |3], |], [5[] and the experimental |6|, [7J communities. Predictions for many of them already 
exist in the literature || |^, 10, 11, |1^, T3, H], and several new experimental investigations are under way 



or planned. In particular, initial data on the process Kl — > 7i°~fe + e~ || have already been presented 
[0], together with further data on Kl — > 7r°77. Confident that the objects of our calculation, the decays 
K + — ► 7r + 7/ + /~, are accessible to experiment, our goal is to provide information on their rate and the 
corresponding decay distributions. 

Analogously to what has been studied in the case of Kl — > 7r°77, the reaction K + — > 7r + 77 takes 
place predominantly through loop diagrams with pions in the loop. In the former process, the decay 
distribution is quite distinctive and the rate is predicted without any free parameters at one-loop order. 
While the distribution agrees well with experiment, the theoretical rate appears too small by more than 
a factor of 2. Because of this, several authors have gone beyond the straightforward one-loop (order 
E 4 ) chiral calculation. Adding a series of higher order effects in a quasi-dispersive framework one has a 
surprising success at increasing the rate without modifying the decay distribution greatly |10], ||. The 



process K + — > 7r + 77 has been studied in a similar way [|T^, 13]. The physics which determines the above 
decays is also involved in the reaction we consider. The experimental study of the lepton-photon modes 
can achieve independent confirmation of the dynamics that drives the whole complex of decay modes. 

As was done in previous works, we separately calculate the one-loop results within ChPT both at 
0(E 4 ) and (9(E 6 ). The first one gives a prediction for the rate and the variation of the amplitude 
depending on the invariant mass of the two leptons, kf, which is carried by an off-shell photon. An 
additional parameter, in the form of a local counterterm, also enters the calculation. In the second case, 



following Ref. |I2" |, we take into account the higher order behavior in the experimental K + — > 3ti decay 
rate. 

At (9(E 6 ) the higher order effects in the K + — > 3tt vertex are extracted from a quadratic fit to the 



amplitude. According to the results of Refs. [O, Oj, we shall not be concerned with vector meson 



corrections, likely too small to be significant, given the uncertainties in the several parameters involved 
in this computation. 
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This paper is organized as follows: In Sec. II we fix our notation and define the quantities used in the 
rest of this paper by summarizing some established results for the decay K + — > 71+77. This provides a 
starting point for our calculation, taking a photon off shell and going through the process K + — > 7r + 77* 
— > 7r + 7/ + /~. In Sec. Ill we describe the 0(E A ) calculation, which we extend to 0(E 6 ) in Sec. IV, taking 
into account the unitarity corrections at one loop. Finally, we recapitulate our conclusions in Sec. V. All 
the relevant expressions for the integrals used in this paper are shown in the Appendix. 

II. K — > 7T77 amplitudes 

Let us first review some previously known results for K — > 7777, and establish our notation for the 
following sections. We define the general amplitude for K — > 7T77 as given by 

M[K{p K ) -> 7r(p) 7 (A; 1 , e l ) 1 (k 2 , e 2 )] = e lfl e 2l/ M^(p K , k u k 2 ) (1) 
where ei,e 2 are the photon polarizations, and M*" has four invariant amplitudes: 
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The physical region in the adimensional variables y and z is given by 
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Here k\ and k 2 are the momenta of the off-shell and on-shell photons, respectively, with the off-shell photon 
materializing into the lepton pair. Note that the invariant amplitudes A(z, y), B(z,y) and Ci(z,y) have 
to be symmetric under the interchange of ki and k 2 as required by Bose symmetry, while C 2 (z,y) is 
antisymmetric. 

Using the definitions the double differential rate for unpolarized photons is given by 
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In the limit where CP is conserved, the amplitudes A and B contribute to K 2 — > 7r°77 whereas K\ 
involves the other two amplitudes C% and C 2 . All four amplitudes contribute to K + — > 7r + 77. Only A and 



(6) 



7T°77 
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C\ are non- vanishing to lowest non-trivial order, 0(E A ), in ChPT. As argued in |15|, [12|, the antisymmetric 
character of the Ci{z,y) amplitude under the interchange of k\ and k 2 means effectively that while its 
leading contribution is 0(E 6 ), this can only come from a finite loop calculation because the leading 
counterterms for the C2 amplitude are 0(E 8 ). Moreover, this loop contribution is helicity suppressed 
compared to the B term. This antisymmetric (9(E 6 ) loop contribution might be smaller than the local 
0(E S ) contribution. 



III. 0(E 4 ) calculation 



First let us provide the straightforward 0(E A ) calculation of Ai(K + — > tt + ~/1 + 1~) within ChPT. This is 
the generalization to k\ 7^ of the original chiral calculation of the authors of ]IE , [17] , and it includes all 
the k\jm\ and q = k\/m? K variations of the amplitudes at this order in the energy expansion. There can 
be further k\/(l GeV) 2 corrections which correspond to (9(E 6 ) and higher. The easiest technique for this 
calculation uses the basis where the kaon and pion fields are transformed so that the propagators have 
no off-diagonal terms, as described in Refs. |TE|, |T?J. Some of the relevant diagrams are shown in Fig. 1. 





a) 



b) 




c) d) 

Figure 1: Some diagrams relevant to the process K + — > n +r yl + l~ at 0{E A ) and 0(E 6 ). Either photon 
may also be radiated from the incoming K + or the outcoming 7r + . The lepton pair must be attached to 
one of the photons, and the on-shell photon may be radiated from one of these leptons. 

Analogously to the leading AI = 1/2 0(E A ) A(z,y) and C\(z,y) amplitudes for K + — > 71+77 which 
have been computed in ]T7 , we can write an expression for the amplitude for K + — > 7r + 77*, 



(z) = 2 ^ e ™ q) {{z + 1 - rl)[l + 21 (ml)] + (z + r\ - 1)[1 + 2I(m\)\ - c(z - q)} , (7) 
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where Gs is the effective weak coupling constant determined from K — > -kit decays at 0(E 2 

G F 
= 5.1, 



where is the Cabibbo-Kobayashi-Maskawa matrix [18], and 
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The above functions are related to those presented in Ref. 
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which are shown below: 
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The above results agree with the results obtained in [12] in the k\ — > limit. 

In Eq. (R]) the pion loop contribution largely dominates over the kaon loop part. The loop results are 



finite, but ChPT allows an (9(E ) scale independent local contribution that may be parametrized as [JT9 
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a = ^^[3(L 9 + Lio) + #14 - #15 - 2JV 18 ] 



(17) 



or, using the notation of [17, 20 



32tt' 



[12(1,9 + L 10 ) -wi- 2w 2 - 2w 4 ], 



where c is a quantity of 0(1). The L 9 and L w are the local C(E 4 ) strong couplings and #14, #15 and 
#8 (or Wi, W2 and W4) are 0(E 4 ) weak couplings, still not completely fixed by the phenomenology, and 
which can be only computed in a model dependent way pi| . The weak deformation model (WDM) |20[ 
predicts c = 0, while naive factorization in the factorization model (FM) |22|, ^TJ gives c = —2.3. In these 
models, because of the cancellation in the vector meson contribution in c, the role of axial mesons could 
be relevant [pl|l . 

Thirty-one events for the process K + — > 71+77 have been observed at BNL (E787) ||23|| , with the partial 
branching ratio BR(K+ 71+77, 100 MeV/c < P n + < 180 MeV/c) = (6.0±1.5{stat}±0.7{syst» x 10~ 7 . 
This has been extrapolated with the help of ChPT, performing a maximum likelihood fit of c to the 
spectrum. The results of the fit to the data support the inclusion of the unitarity corrections, giving as 
the best fit c = 1.8 ± 0.6 and BR(K + -> 7^77) = (1.1 ± 0.3 ±0.1) x 10" 6 , as also reported in the 
Review of Particle Physics |24| . 

The C(E 4 ) contribution to the Ci(z,y) amplitude is 



Ci(z) 
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7T 
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2 + rt 



z -r 7T + z'jv 



z — r r 



m K 
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where r v = m v /rriK and T 7 _ 



7T" 



Witten functional [f25| 



77 



0. This amplitude is generated by the Wess-Zumino- 



^7r", 77) — > 77* through the sequence K + — > 7r + (7r u , 77) — > 7r + 77*. This contribution 
amounts to less than 10% in the total width. 

The (9(E 4 ) results can be expressed as total branching ratios. They are summarized in Table I for 
three values of c, given respectively by the weak deformation model, the factorization model, and the fit 
to BR(i^ + — > 7r + 77) mentioned above. 



Table 1: Results for BR(K+ -»■ ir+jl+l-) at C(E 4 
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7/x + /z ) 
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c = (WDM) 


8.6 
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3.6 
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c = -2.3 (FM) 


5.7 
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10 


-9 


3.9 


X 


10 


-11 



The decay distributions in z and y provide more detailed information. We present them in Figs. 2-5. 



IV. 0(E 6 ) calculation 



In this section we extend this calculation along the lines proposed by the authors of Refs. 0, [T(| for 



decays and D'Ambrosio and Portoles for K + decays |12| . The former provided a plausible solution to the 
problem raised by the experimental rate not agreeing with the 0(E A ) calculation when both photons are 
on shell. We have to add a new ingredient that involves known physics that surfaces at the next order in 
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Figure 2: The differential branching ratio 
dBK(K + — > 7r + 7e + e _ )/dz to order E 4 is plot- 
ted vs z for c = 1.8 (solid line), c = (dashed 
line) and c = -2.3 (dotted line). 



Figure 3: The differential branching ratio 
dBR(i^ + — > 7r + 7e + e~)/dy to order E 4 is plot- 
ted vs y for c = 1.8 (solid line), c = (dashed 
line) and c = -2.3 (dotted line). 



dBR(K + ^TT + 7yti + / Li")/dz 




z 



Figure 4: The differential branching ratio 
dBR(fT + — > 7r + 7// + //~)/dz to order E 4 is plot- 
ted vs z for c = 1.8 (solid line), c = (dashed 
line) and c = -2.3 (dotted line). 



dBR(K + ^77 + 7Ai + /i,")/dy 




y 



Figure 5: The differential branching ratio 
dBK(K + — > 7r + 7/i + /!~)/dy to order E 4 is plot- 
ted vs y for c = 1.8 (solid line), c = (dashed 
line) and c = -2.3 (dotted line). 
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the energy expansion, i.e. the known quadratic energy variation of the K 
from higher order terms in the weak nonleptonic Lagrangian |26|, [K| 0. 



27 



3ir amplitude, which occurs 
While the full one-loop 



structure of this is known 



it involves complicated nonanalytic functions and we approximate 



the result at 0(E ) by an analytic polynomial which provides a good description of the data throughout 



the physical region [p0| , ^q] . Expanding in powers of the Dalitz plot variables, 



2(Ci + Ca) {y 2 + ^) " (6 + Ca - Ca) " ^) • 



(20) 



Here the subscripts 1 and 3 refer to AI = 1/2,3/2 transitions respectively, and the coefficients in Eq. 
fl20"D have been fitted to the dat a p8| . We omit the AI = 3/2 couplings £3 and £3, £ 3 because of their big 
errors shown in the fits in Ref. ||28||. The Dalitz plot variables are commonly defined as 

X = Y = (21) 

7T 7T 

with Si = (pk — Pi) 2 for i — 1, 2, 3, so = (s\ + s 2 + «3)/3, and the subscript 3 indicates the odd pion (n° 
for K L decays and n~ for K + decays). 

In principle one can add the ingredients to the amplitudes and perform a dispersive calculation of the 
total transition matrix element. In practice it is simpler to convert the problem into an effective field 
theory and and do a Feynman-diagram calculation which will yield the same result. We follow this latter 
procedure. 

The Feynman diagrams are the same as shown in Fig. 1, although the vertices are modified by the 
presence of 0(E 4 ) terms in the energy expansion. Not only does the direct K — > 3tt vertex change to the 
form given in Eq. (p0[) , but also the weak vertices with one and two photons have a related change. The 
easiest way to determine these is to write a gauge invariant effective Lagrangian with coefficients adjusted 
to reproduce Eq. (^). One also has to add diagrams with one or two photons radiating from the incoming 
K + or the outcoming 7r + , or with one photon radiating from one of the outcoming leptons. This causes 
infrared divergences, which are to be treated in the usual way, as part of a general calculation of radiative 
corrections to the process K + — > 7r + / + / _ . In practice, with an appropriate set of experimental cuts on 
the phase space parameters, it is possible to restrict the outcome to a measurable non-bremsstrahlung 
contribution only H31J| . Below we shall give an example of these cuts and a prediction for the experimental 
result once they are implemented. 

The resulting calculation follows the same steps as described in Sec. Ill, but is more involved and is 



not easy to present in a simple form. We have checked that our result reduces to that of Ref. [0 in the 
limit of on-shell photons. Remembering the definitions 

f n = , r„ , z q (22) 

m K m K m L K m A K 

the unitarity one-loop corrections yield the following: 
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where 
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1 



3^(2Ci - 6)^ [1 + 6(r 2 - *) + 9(r 2 - zf\ j [1 + 21 (ml)] 
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(4Ci + 6) - {2[1 - 2(xi + rr 2 )]/i(^ 2 ) + 0:1/1(22) 
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The integrals used in the above formulas are defined here and given explicitly in the Appendix: 
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The above formulas lead to the total branching ratios shown in Table II, in full analogy with the 
results of Sec. III. The numerical results are obtained for the mass scale fi = m p and setting all the 
counterterms to 



12]. The corresponding decay distributions are plotted in Figs. 6-9. 



Table 2: Results for BR(K + -> 71+7/+/-) at C(E 6 ). 
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7/i+/i ) 


c = 1.8 (fit) 


1.7 
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7.0 
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c = (WDM) 
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7.3 
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10 
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c = -2.3 (FM) 


9.2 


X 


10- 9 


8.5 
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10 


-11 



The uncertainty in the theoretical prediction is dominated by the unknown 0(E 4 ) counterterm gener- 
ated amplitude c in Eq. (0). In Figs. 10 and 11 we plot BR(fT + — > n + je + e~) and BR(K + — > n + jfi + fi~) 
as a function of c, both with and without the 0(E 6 ) corrections just computed. 

If we implement the experimental cuts currently used at BNL to extract the non-bremsstrahlung 
contribution to K + — ► 7r + 7e + e~ pi| , 

m e+e - > 150 MeV, E 1 > 30 MeV, m e+7 ,m e - 7 > 30 MeV, (34) 

the resulting theoretical branching ratios are reduced by more than an order of magnitude. They are 
presented in Table III. Preliminary experimental data are not conclusive at the present stage. In Figs. 
12 and 13 we plot the differential branching ratios up to (9(E 6 ), taking into account the above cuts in 
the phase space integration. 
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Figure 6: The differential branching ratio 
dBR(A" + -> 7T + 7e+e-)/dz to order E 6 is plot- 
ted vs z for c = 1.8 (solid line), c = (dashed 
line) and c = -2.3 (dotted line). 

dBR(K + ^7T + 7 / u + y u")/dz 




z 



Figure 8: The differential branching ratio 
dBR(A+ -> 7r + 7/i + /i-)/dz to order E 6 is plot- 
ted vs z for c = 1.8 (solid line), c = (dashed 
line) and c = -2.3 (dotted line). 



dBR(K + ^7T + 7e + e")/dy 




-0.50 -0.35 0.00 0.25 0.50 

y 



Figure 7: The differential branching ratio 
dBR(A" + -> 7r+7e+e-)/dy to order E 6 is plot- 
ted vs y for c = 1.8 (solid line), c = (dashed 
line) and c = -2.3 (dotted line). 



dBR(K + ^7T + 7Ai"V~)/dy 




y 



Figure 9: The differential branching ratio 
dBR(A+ -> 7r+7 / u+/i-)/dy to order E 6 is plot- 
ted vs y for c = 1.8 (solid line), c = (dashed 
line) and c = -2.3 (dotted line). 



Table 3: Results for BR(K + -> 7r + 7e+e-) at C(E 6 ) with the cuts defined by @. 



BR(A+ -» 7r+7e+e-) 



c = 1.8 (fit) 5.5 x 10- 

c = (WDM) 4.8 x 10- 

c = -2.3 (FM) 4.7 x 1Q~ 
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BR(K + ^7T + 7e + e") 




-4-2 2 

C 



Figure 10: The branching ratio BR(i^ + — ► 
7T + 7e + e~) is plotted vs c at 0(E A ) (dashed 
line) and up to (9(E 6 ) (solid line). 



BR(K + ^tt + 7aiV") 




6.010 -11 — 
4.010 -11 — 



2.010 -11 — 



-4 -2 2 4 

C 

Figure 11: The branching ratio BR(f^ + — ► 
7r + 7/i + yU _ ) is plotted vs c at (9(E 4 ) (dashed 
line) and up to 0(E 6 ) (solid line). 



dBR(K + ^7T + 7e + e )/dz 



dBR(K + ^7T + 7e + e )/dy 




5.010 -10 




Figure 12: The differential branching ratio 
dBR(K + -> 7r+7e + e")/dz to order E 6 is plot- 
ted vs z for c = 1.8 (solid line), c = (dashed 
line) and c = -2.3 (dotted line) with the cuts 
defined by the inequalities 



Figure 13: Same as in Fig. 12 for the differen- 
tial branching ratio dBR(K + — > 7T + 7e + e - )/dy 
vs y. 



11 



V. Conclusions 



We have computed the unitarity corrections at one loop in ChPT for the processes K + — > 7r + 7e + e~ and 
K + — > 7r + 7yU + yU _ , allowing us to present predictions for their total and differential branching ratios. 

As expected, the muonic rate is significantly smaller than in the corresponding electronic mode, 
analogously to the cases studied in Refs. and j|. This is of course due to the more limited phase 
space, as well as the fact that the photon propagator is further off shell in the muonic case. We again 
see that the more complete calculation presented above leads to a conspicuous enhancement over the 
purely order E 4 calculation presented first. The vector meson diagrams here are not expected to add 
a significant amount to the overall rates and have been omitted altogether. Their inclusion provided a 
consistent contribution to the Kl decays, leading to dramatic increases in the 0(E 6 ) results. Therefore, 
in the computation presented in this paper the enhancement is expected to be somewhat smaller. 

The results for the differential branching ratios follow a pattern recognizable in all the previous 
calculations for radiative rare kaon decays: A large peak is visible above the two-pion threshold in the 
z variable, with a tail extended to low z, and a slightly asymmetrical and structureless distribution in 
the y variable. Experimentally, the abundance of information supplied by the former plots makes them 
a preferred option (see for example [0]). This is true also if we reduce the phase space of integration 
performing experimental cuts. 

We also note that comparing Figs. 10 and 11, one sees that the dependence on c is markedly higher 
when electrons, instead of muons, are present in the final state. This has important consequences if one 
tries to extract the value of c from the data. 
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Appendix: Relevant integrals 

In this appendix we list the explicit expressions for the integrals used in the calculation of Sec. IV. We 
follow the notation of that section. For s < Ami and k\ < Ami we have 



s-k( 



kl 



F(s) - F(kl] 



Ami ~ k\ arctan , 



kl 



\ A^kJ + V Aml ~ S arCtan V j ' 
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As (Ami ~ s ) 
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